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The increased flexibility of long endurance aircraft having high aspect ratio wings necessitates attention 
to gust response and perhaps the incorporation of gust load alleviation. The design of civil transport aircraft 
with a strut or truss-braced high aspect ratio wing furthermore requires gust response analysis in the tran- 
sonic cruise range. This requirement motivates the use of high fidelity nonlinear computational fluid dynamics 
(CFD) for gust response analysis. This paper presents the development of a CFD based gust model for the 
truss braced wing aircraft. A sharp-edged gust provides the gust system identification. The result of the sys- 
tem identification is several thousand time steps of instantaneous pressure coefficients over the entire vehicle. 
This data is filtered and downsampled to provide the snapshot data set from which a reduced order model is 
developed. A stochastic singular value decomposition algorithm is used to obtain a proper orthogonal decom- 
position (POD). The POD model is combined with a convolution integral to predict the time varying pressure 
coefficient distribution due to a novel gust profile. Finally the unsteady surface pressure response of the truss 
braced wing vehicle to a one-minus-cosine gust, simulated using the reduced order model, is compared with 
the full CFD. 
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Filtered pressure coefficient 

Predicted value of pressure coefficient 

Unfiltered pressure coefficient 

Gust function factor based on aircraft weight ratios 

Generalized force due to gust 

Gust length, ft 

Free stream Mach number 

Number of structure modes 

Number of surface points 
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Number of time steps 

Number of proper orthogonal decomposition modes 
Free stream dynamic pressure, psi 
Integration matrix, in 2 
Time, sec 

Non-dimensional time 
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Free stream velocity, ft / sec 

z-direction gust velocity magnitude, ft / sec 

z-direction one-minus-cosine gust intensity, ft / sec 

z-direction sharp-edged gust intensity, ft / sec 

z-direction reference gust intensity, ft / sec 

Wave number 

Vector of time coefficients in proper orthogonal decomposition analysis 
Predicted value of /3 
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<5 Structure mode shape matrix 

<P Singular vectors of a proper orthogonal decomposition 

E Vector of training values of dynamic component of surface 

nodal pressure coefficients 
E Predicted value of E 


I. Introduction 

Aerodynamic efficiency increase and drag reduction are key NASA Subsonic Fixed Wing Program goals. 1 A 
component of that research program investigates the truss braced wing (TBW) configuration. Multidisciplinary de- 
sign optimization (MDO) studies of truss braced wing airplanes suggest that optimal designs can have very flexible 
wings. 2 It is well known that vehicles with long flexible wings can have aeroelastic issues related to flutter, gust loads, 
maneuver loads, limit cycle oscillation, ride quality and buckling. 3 The TBW may also have stability issues with low 
sweep angles and very low structural frequencies that can couple with aircraft rigid body modes and a flight control 
system. 

Other aircraft development efforts aimed at high endurance have resulted in very flexible wings. Experience in 
aircraft design programs intended to improve aerodynamic efficiency such as the HiLDA (High Lift-to-Drag Active) 
wing, 4 high altitude long endurance (HALE), 5 and the Aeroenvironment Helios crash investigation 6 indicate that gust 
response requires more thorough analysis and validation using nonlinear multidisciplinary aeroservoelastic codes. For 
these reasons it is likely that many future projects will necessitate a new level of analysis not seen in current production 
aircraft design practices. Final designs will likely include closed-loop flutter suppression and gust load alleviation. 7,8 

Gust analyses have been a standard part of vehicle loads analysis for many decades. Analysis to date indicates 
that gust loads and closed-loop gust response of the TBW may be a critical design factor. Production vehicle design 
practice uses gust analysis with linear aerodynamics. Very sophisticated, but fully linear gust models, have been 
developed for both time and frequency domain analyses. The Laplace transform of an arbitrary gust wave form has 
facilitated the development of frequency domain models. Alternately, a reduced order model (ROM) developed using 
time history data of an appropriate parameter set for an aeroelastic model can be imported into a linearized state- space 
model for closed-loop analysis. However, the reduced order gust model and the aerodynamic response data needed 
to construct that model must be generated first. Time domain gust analysis has historically been performed using a 
panel code in which the introduction of a perturbation velocity as a local angle of attack increment is relatively straight 
forward. A linear panel code is acceptable if the flow field and unsteady aerodynamic response are entirely linear. 
Linear gust aerodynamics may or may not be adequate as more flexible vehicles are designed to fly in the transonic 
flight regime. If the steady state or unsteady flow are nonlinear, a higher order CFD simulation of the gust response 
will be required. 

For this reason there has been an interest in developing gust modeling capability in high fidelity nonlinear CFD 
codes. A technique called the field velocity method (FVM) has been developed and is widely being implemented for 
the simulation of a gust within CFD codes. 9-12 A previous publication has discussed the development of this method 
in the NASA CFD code FUN3D. 13 The present paper uses the gust simulation capability in FUN3D to perform gust 
analysis and develop a gust model for the truss braced wing aircraft. The development of a CFD based gust model also 
provides an opportunity to explore a new approach to reduced order gust modeling. A method of creating the ROM of 
a gust is discussed using proper orthogonal decomposition (POD) combined with a convolution of the gust response. 

At the same time an accurate gust model based on CFD will provide value added to the truss braced wing (TBW) 
aircraft design effort. Development of that gust model is the focus of the present paper. The TBW version to be used 
in this paper is the Subsonic Ultra Green Aircraft Research (SUGAR) High aircraft being developed by the Boeing 
company. 14 This vehicle is being designed under a NASA Fixed Wing N+3 NRA contract. This effort is tasked with 
developing an aircraft that might carry passengers in the next 30 year time frame. By using truss structures the aircraft 
can be designed with a thinner and lighter inboard wing, thus reducing weight and making possible the use of a higher 
aspect ratio wing than is used by current civil transports. The increased lift to drag ratio that is anticipated contributes 
to a reduced fuel usage. The intent in this paper is to develop a reduced order gust model that can be used in a closed- 
loop aeroservoelastic design and analysis of a gust load alleviation control law. The first sections of this paper will 
briefly discuss the FUN3D code, the CFD model, the field velocity method of modeling gusts, and the approach to 
creating the reduced order model. The reduced order model results will be compared with a direct FUN3D simulation 
of a one-minus-cosine gust. 
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II. FUN3D Solver 


The Navier-Stokes code FUN3D (fully unstructured Navier-Stokes three-dimensional) is a finite-volume unstruc- 
tured CFD code for either compressible or incompressible flows. 15,16 Flow variables are stored at the vertices of the 
grid. FUN3D can solve the discrete compressible Euler or Reynolds-averaged Navier-Stokes (RANS) flow equations 
either tightly or loosely coupled with a turbulence model on mixed element grids, including tetrahedra, prisms, pyra- 
mids and hexahedra. The inviscid fluxes are computed using an approximate Riemann solver based on the values 
on either side of a cell interface. Cell interface values are obtained by a least squares extrapolation using gradients 
computed at the vertices. 

The present study uses the RANS solution capability of FUN3D. The simulations include turbulence provided 
by loosely coupling the Spalart-Allmaras turbulence model. 17 In the present study the low dissipation flux splitting 
scheme for the inviscid flux construction and the Venkatakrishnan limiter are used. The inviscid flux and viscous terms 
are computed with second order accuracy. The unsteady solutions at each time step are updated with a second order 
backwards Euler time differencing scheme. Steady state and subiterative solutions are accelerated to convergence by 
the use of local time stepping. 15 Domain decomposition is used to enable distributed parallel computing. 


III. The Field Velocity Method and Discrete Gust Model 


The FVM physically introduces gusts into a CFD code by utilizing grid velocity. 9-12 Normally grid velocity would 
be associated with physically moving the grid. However, it is possible to introduce an arbitrary perturbation velocity in 
a stationary grid by prescribing the grid velocity at either all or some of the field grid points without actually moving 
the grid. For instance, the gust profile can be defined by a streamwise variation in the z-direction perturbation velocity 
such as 

Wg(x,t) =f z (t-x/Uoo) ■ (1) 

This method has been used for relatively simple to moderately complex configurations, for example Singh and 
Baeder, 11 Zaide and Raveh, 18 Raveh, 19,20 and Wang et al. 21 It was recently used to simulate the aeroelastic response 
of a launch vehicle to a sequence of one-minus-cosine gusts. 13 

In this study a one-minus-cosine and a sharp-edged gust profile will be used. In the present case, having only a 
z-component, the one-minus-cosine gust velocity is 


w g (x,t) 


_ 1 
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1 — COS 


( xjt-x/uj Y 

V h/u- /. 


for tUoo — 2 L g < x < tUoo 


( 2 ) 


where L g is the gust length and w omc is the one-minus-cosine gust intensity. In order to obtain the reduced order (ROM) 
gust model, a time-accurate unsteady solution is computed for the response to a sharp-edged gust. The sharp-edged 
gust velocity is given by 


w g {x,t) = wseg for t > x/Uoo and w g (x,t) = 0 for t<x/Uoo . 


( 3 ) 


IV. Development of a POD/Convolution Based Reduced Order Gust Model 

Rather than computing a fully time-accurate CFD solution for each gust profile, it would be of interest to create 
a reduced order model that is capable of predicting the time history response of the distributed surface pressures to 
a gust input. Such a model can be constructed using the convolution of a sharp-edged gust, combined with a model 
constructed from the proper orthogonal decomposition of the covariance matrix of the sharp-edged gust unsteady 
pressure coefficients. In this model, the POD eigenvectors provide the spatial distribution of the principle components 
of the unsteady pressure response. Combining the POD eigenvectors with a convolution integral provide a method of 
predicting the unsteady surface pressure responses to an arbitrary gust. 

A sharp-edged gust (SEG) velocity profile is used to generate the pressure response training data. For an array of 
surface forces that are continuous in time, the aeroelastic generalized forces due to a gust can be written as 

GsEc(t) = [8} T f(t) (4) 

where / are the aerodynamic nodal forces derived from the pressures at CFD surface grid points and [8] are the 
structural mode shapes. The time history of generalized forces due to an arbitrary gust input can be constructed by 
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convolution, namely 


G g (t) = w g (T)G sEG (t -z)dz . (5) 

As Raveh does, it is possible to compute the generalized force for any arbitrary gust using the generalized force Gseg 
due to the sharp-edged gust. 20 On the other hand, it is useful to compute the convolution using a POD model of the 
unsteady pressures. By using a POD, the nodal pressure coefficients due to the sharp-edged gust can be decomposed 
into the spatial component of the POD <f> and the time varying component /3(f). Due to an arbitrary gust, the nodal 
forces are 

f(t) = q°°[s][®]P(t) (6) 

(the hat indicates the predicted value) and the generalized forces are 

G g (t) = [P}J‘w g (t)l3 SEG (t-t)dt=[P}p(t) (7) 

where [P] = ^[S] 7 ^] [<£>]. 

The discrete time formulation of the POD/convolution requires additional discussion. A proper orthogonal decom- 
position of the surface pressures for all time steps of an entire CFD simulation is out of the realm of possibility except 
for a very coarse mesh on a simple geometry. The present aircraft mesh has nearly 500,000 surface nodes, which by 
today’s standards is not an excessively refined grid. At the present refinement, it has been found necessary to use a 
subset, or snapshots of the surface pressure time history. The POD eigenvectors are computed by using the matrix of 
pressure coefficient data at N ss snapshots 

[X] r =[S! ••• E, ••• (8) 


The columns of this matrix are composed of the surface node pressure coefficients at a given snapshot. The snapshot 
/ of the training data is 


c P\,l 




L C PN s ,l\ 


(9) 


The relationship between the array of pressure coefficients and the POD model at time step / is given by £/ = [<t>] j 3j . 
The matrix <f> is the collection of singular vectors of the matrix [X ] T . The size of the matrix d> is reduced to form 
a truncated POD-basis <F r with dimension N s x Npod by retaining only the eigenvectors associated with the largest 
eigenvalues. Although it is possible to find the eigenvectors of the covariance matrix, in practice, especially for large 
systems, computing the eigenvectors from a singular value decomposition (SVD) is more efficient and robust. The 
rank-r approximation of the matrix X T is obtained as an approximate solution to 


X T « 


( 10 ) 


by an SVD. The orthonormal vectors y /\ . 1 // 2 , •••, y/ r are the approximate left singular vectors and the orthonormal 
vectors , $ 2 , are the approximate right singular vectors, while G\ > 02 >,..., cr r are the singular values of X T . 
The matrix <f> r contains the leading r right singular vectors of the matrix X T . A very efficient stochastic algorithm 
for computing the singular value decomposition of large systems has been developed in Liberty et al. 22 and Rokhlin 
et al. 23 This algorithm will be used and is discussed further in the next section. Finally, if the POD modes are unit 
normalized so that <F r <F = /, the time varying component of the POD model of the training data can be obtained for 
time step / by p SEGl i = [®r] T Zi- 

The N ss snapshots of training data are a subset constructed by taking CFD surface pressure coefficients at constant 
intervals over the entire sharp-edged gust simulation. The downsampling of a data set to obtain snapshots requires an 
anti-aliasing filter so as not to corrupt the data and potentially destablize a numerical scheme with unwanted higher 
frequencies. Filtering does not appear to be the standard practice in the literature when creating POD models from 
CFD pressure data. It would seem to be something of growing importance, however, as more CFD models include 
detached eddy simulation (DES) or large eddy simulation (LES) turbulence models and as larger grid sizes require 
data downsampling. The filtering used here is by Lele. 24 The filtered surface pressure data is computed using the 
sixth- order scheme 

D C B 

bCp,i— 2 MCp,i— 1 + Cp i + CLCp i-\-\ ~\~ bCp i+2 — Aff + — (Cp^-i-3 + Cpj— 3 ) H - ^ (Cp,i+2 H“ Cp,i—l) H - ^ (Cp,i+ 1 Cp,i— 1 ) (1 1) 
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The c p j are the unfiltered pressure coefficients while the c p j are the filtered values for i : 1 N ss . The coefficient 
a = 0.4999 was chosen to give the best behavior and b — (3 — 2a)/ 10, A = (2 + 3<z)/4, B = (6 + 7<z)/8, C = (6+a)/20, 
D = {2 -3a) /40. 

The convolution is integrated with a discrete method of the form 

/-t 

fc=l 

where the coefficients c * are chosen to give second order accuracy. Surface pressures can be reconstructed for an 
arbitrary gust by 

S/ = [<*V]A • (13) 

V. Stochastic Algorithm for Computing Singular Value Decomposition 

As the size of CFD models becomes larger, there is a growing need to implement memory and computing efficient 
algorithms for performing the proper orthogonal decomposition. There are several methods available to perform the 
factorization in equation 10, however, a very efficient and robust method is the stochastic algorithm of Liberty et al. 22 
and Rokhlin et al. 23 and Halko et al. 25 The source code developed by those authors, and available from the home page 
of one of the authors, is used in the present computations. The brief description of the algorithm is given here that 
follows Martins son. 26 Let A = X T . Sample matrices can be computed 

Y =A£l and Z = A T 0 (14) 

where Q and 0 are random matrices. Orthonormal matrices Q and W can be found so that 

Y = QQ T Y and Z = WW T Z . (15) 

Solve the linear systems 

Q t Y = T(W t &) and W r Z = T r (Q r 0) (16) 

for T. Factor T such that 

T = UIV T . (17) 

The desired orthonormal matrices are found by setting 

U = QU and V = WV . (18) 

The matrix V contains the desired eigenvectors. This method has been successfully used for extremely large systems. 
Previous studies by the authors have shown it to be very efficient from a memory and CPU usage standpoint compared 
to standard methods of singular value decomposition. 


VI. Results 

The TB W outer mold line used in the present analysis is shown in Figure 1 . The CFD volume mesh was generated 
using VGRID 4.00. 27 It is a half-span tetrahedral mesh with 17.3 million nodes with mirror boundary conditions at the 
vehicle symmetry plane. The surface is modeled with 478,000 nodes. The normal spacing at the wall is approximately 
at y+ = 1 everywhere. Figure 2 shows the surface mesh clustering in the wing and engine nacelle regions. This paper 
will investigate the development of a gust model that can be used in a future closed-loop aeroervoelastic analysis of 
the truss braced wing aircraft. The focus of the present analysis is a flight condition at which Boeing analyses have 
resulted in the largest gust induced wing bending moment. That gust is a one-minus-cosine gust with gust length of 
1 10 ft and gust intensity of 31.4 ft/sec. The Boeing aeroelastic analysis was done with MSC.Nastran™ using doublet 
lattice aerodynamics and a linear structural model. The condition at which that analysis was done is Mach 0.70 at 
an altitude of 19k ft and zero degrees angle of attack. The present rigid vehicle gust analysis uses the same flight 
condition. 

The gust simulations are initialized with a steady state solution. Steady state solutions were obtained after 30,000 
iterations using Courant-Friedrichs-Lewy (CFL) number of 10 for the flow and 6 for the turbulence equations. The 
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steady state lift coefficient varied less than 10 -5 over the last 5000 iterations. Likewise, the moment coefficient varied 
less than 10 -4 over the last 5000 iterations. Both the lift and moment coefficients had jitter that continued to the end, 
but that was also extremely small in comparison to absolute values of the coefficients. This jitter had been attenuated 
somewhat by the use of the Venkatakrishnan limiter but was still present. That jitter also shows up minutely in the 
unsteady pressure coefficient time histories. Returning to the steady state solution, the surface pressure distribution 
on the aircraft is shown in Figure 3(a). Note that the model includes engine nacelle with intake and both fan and 
core exhaust planes. Inlet and exhaust conditions are set in the present computations to produce reasonably realistic 
flow conditions around the engine. Contours of stream wise velocity around the engine are shown in Figure 3(b). The 
exhaust flow from the engine core is evident from the colored contours in that figure. 

For the unsteady simulations a dimensional time step of 3.2 x 10 -4 sec was used. Twenty sub-iterations per time 
step with a CFL number of 15 for the flow and 8 for the turbulence model were used. Pressure coefficient data was 
written at each time step. The complete time history of data included 3000 time steps. This many data samples with a 
surface of nearly 500,000 nodes proved to be too large a data set to reasonably compute the S VD on a single computing 
node. This computing limitation required that the pressure coefficient data be downsampled to something less than 
500 snapshots. For this reason, the surface pressure data was repeatedly filtered and downsampled by two with the 
method described earlier, to a subset of every eighth point in time. This left a data sample with 376 snapshots. An 
alternate downsampled data set was also extracted that was unfiltered. The POD analysis was next performed on the 
two downsampled data sets using the stochastic S VD software by Liberty et al. 22 and Rokhlin et al. 23 Models with 
10, 20, 30 and 40 POD modes were developed. It should be noted that a standard method of SVD was also used with 
which to compare. The stochastic algorithm easily out performed the standard SVD algorithm in execution time. The 
solutions by the two methods were nearly identical. 

A global error analysis was performed by comparing the predicted surface pressure time histories with the original 
surface pressure time histories due to the sharp-edged gust. A plot of the global error from the filtered and unfiltered 
data sets as a function of the number of POD modes is shown in Figure 4. It can be noted that the addition of modes 
significantly reduces the error, at least in the case of the filtered data. The unfiltered data error reduces initially, but 
clearly the return diminishes as the number of modes increases. 

A comparison of the time histories of the ROM predicted and FUN3D simulated surface pressures at a collection 
of points will be shown for the one-minus-cosine gust. Figure 5 shows the point stations at which time histories will 
be presented. Several stations are chosen on the upper wing surface, several on the upper surface of the strut, one on 
the engine nacelle and one on the horizontal tail. Figures 6(a)-(f) show surface pressure coefficient time histories. 
The values plotted are c' p = ( c p — c Pss )/c p 0 where c Pss is the steady state value. The parameter c po is a normalization 
factor that is used consistently throughout the time histories. Time histories are shown for the FUN3D simulation and 
for the 40 POD mode convolution responding to the one-minus-cosine gust. Station 1 and 6 are for the wing upper 
surface. The c' p time histories at those stations indicate a strong negative pressure dip as the gust passes. Station 2 is on 
the horizontal tail. The time delay of the pressure response at the tail is evident by comparing Figure 6(b) with (a) or 
(f). The strut responses are found in Figure 6(c) and (d). As expected the flow field does not respond as strongly in the 
strut region as it does over the wing. However, strut station 4 in Figure 6(d) clearly is seeing the influence of the under 
side of the wing that induces a positive pressure response. As can be seen in the figure, the FUN3D and the predicted 
time histories compare very well for all the stations shown. Time histories of the lift and moment coefficients given 
by the 40 mode ROM model and the FUN3D simulation for the one-minus-cosine gust are shown in Figures 7-8. 
The values plotted are = (Q _ Q W )/Q 0 where C% ss is the steady state value, C^ Q is a normalization factor and 
t; = (L,M) for lift and moment coefficient. The lift and moment coefficients are computed using equation 7. The 
perturbation due to the gust is modeled very well by the present ROM model. Both the pressure coefficient and the lift 
and moment coefficient time histories show what appears to be drift in the FUN3D solution, as seen toward the end 
of the time trace. Since the FUN3D solution is fully nonlinear, it may be that there is some drift in the solution not 
evident in the linear ROM model. 

It is of interest to know to what degree the gust simulated here is attenuated by the filtering of equation 11. To 
investigate the behavior of the filter, the transfer function of equation 11 can be written, following Lele 24 


2 » 


A + Bcos(w ) + Ccos(2w ) + Dcos(3w) 
1 + 2acos(w) + 2bcos(2w) 


where w = 2 nk/N 


(19) 


and where w is the wave number. The transfer function of the filter (equation 11) for the values used here is shown in 
Figure 9 for the range of wave numbers up to the Nyquist frequency. The wave number based on the one-minus-cosine 
gust length of 1 10 ft is w « 0.053 for the present data set. The gust response fundamental frequency falls well within 
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the range (0 < w < 1.5) in which signals are returned essentially unattenuated. It can also be seen from the figure 
that at the maximum wave number and above, signals are strongly attenuated, as desired. For reassurance one can 
compare the pressure coefficient time histories for the filtered and unfiltered ROM simulated gust. Figure 10 presents 
a comparison of the filtered and unfiltered pressure coefficient time histories at station 6 as reconstructed by the present 
ROM. The filtered and unfiltered responses are nearly identical. In the future, it may be possible to find an anti-aliasing 
filter that is better fitted in general for that purpose, but for the present one-minus-cosine gust, the current filter seems 
to be adequate. 

A Federal Aviation Regulation ( FAR) family of gusts required for certification has the form 

Womc = Wre fF g ( A- ) 1 /6 (20) 

where w re f is specified based on altitude and F g is a factor based on aircraft weight ratios. Having a reduced order 
gust model enables an analysis of the effect of parametric variations of gust parameters on the vehicle response. The 
variation of gust response with gust length L g has been calculated for the truss brace wing aircraft at Mach 0.70 and 
19 k ft. A family of lift and moment coefficient responses are shown in Figures 11(a) and (b). The lines for L g = 110 
ft correspond to the gust giving maximum wing bending moment for the flexible vehicle according to current linear 
analyses. The lift and moment coefficient have been normalized to the response for this gust length. For the rigid 
aircraft, longer gust lengths result in larger vehicle lift and moment responses but also longer time over which loading 
builds. Of course, to get the actual vehicle response and importantly wing bending moment, the effect of vehicle 
flexibility will be required as well. 


VII. Conclusions 

A computational fluid dynamics based gust model has been developed for the Boeing truss braced wing aircraft. 
A sharp-edged gust and a one-minus-cosine gust have been simulated using the computational fluid dynamics code 
FUN3D. The unsteady pressures from the computational fluid dynamics sharp-edged gust simulation are used to create 
a reduced order model from which the surface pressures of any arbitrary gust can be reproduced. A proper orthogonal 
decomposition enables a reduction in the size of the reduced order model. A convolution of the sharp-edged gust 
response enables a prediction of the response to an arbitrary gust. The combined proper orthogonal decomposition 
and convolution reduced order model enables the prediction of the distributed surface pressure time history due to 
an arbitrary gust. As the size of computational fluid dynamics models increases, algorithmic efficiency will become 
increasingly important. A stochastic singular value decomposition algorithm is used to calculate the principle com- 
ponents of the aerodynamic response. This algorithm is found to perform well for the large data set used here. Even 
with the current moderate number of surface nodes, downsampling of the unsteady pressure data was required. A 
digital anti-aliasing filter was used in the filter/downsample step. For the present one-minus-cosine gust, the filtered 
and unfiltered pressure responses were found to be nearly identical. As computational fluid dynamics models increase 
in size (requiring additional downsampling) or as detached eddy or large eddy simulation turbulence models are used, 
it is anticipated that filtering will become an important topic. The reduced order and the FUN3D modeled surface 
pressures and generalized force response to the one-minus-cosine gust are found to agree well. The reduced order gust 
model is used to compute a family of one-minus-cosine gust responses for the rigid truss braced wing aircraft. It is 
anticipated that the surface pressure responses from these simulations will be used to construct a closed-loop reduced 
order model of the flexible truss braced wing aircraft. 
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(a) Wing surface mesh 


(b) Nacelle surface mesh 


Figure 2. Surface mesh. 



(a) Surface pressure coefficient contours 



(b) Streamwise velocity contours 


Figure 3. Steady state contours, Mach 0.70, a = 0 deg. 
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Figure 4. Global error of proper orthogonal decomposition of the sharp-edged gust. 



Figure 5. Stations at which c p time histories are plotted. 
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c 


p 



(a) Station 1, wing upper surface (b) Station 2, horizontal tail 



(c) Station 3, wing strut (d) Station 4, wing strut 



(e) Station 5, engine nacelle (f) Station 6, wing upper surface 

Figure 6. Pressure coefficient time histories due to one-minus-cosine gust (ROM uses equation 13 with 40 modes derived from filtered 
data). 
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t 

Figure 7. Lift coefficient time history due to one-minus-cosine gust (ROM uses equation 7 with 40 modes derived from filtered data). 



t 

Figure 8. Moment coefficient time history due to one-minus-cosine gust (ROM uses equation 7 with 40 modes derived from filtered data). 
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T(w) 



Figure 9. Transfer function versus wave number for filter (equation 19). 



t 


Figure 10. Pressure coefficient time histories due to one-minus-cosine gust, station 6. 
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(a) Lift coefficient 


(b) Moment coefficient 


Figure 11. Variation of lift and moment coefficient time histories as function of gust length. 
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